--AD-A166  189  BAVESIAN  ANALVSIS  OF  INSPECTION  SAMPLING  PROCEDURES 
DISCUSSED  BV  DEHINGCU)  CALIFORNIA  UNIV  BERKELEV 
OPERATIONS  RESEARCH  CENTER  R  E  BARLOU  ET  AL  JAN  86 
UNCLASSIFIED  ORC-86-2  DAAG29-85-K-0288  F/G  12/1 


o> 

00 

<0 

<0 

operations 

5 

RESEARCH 

1 

0 

< 

CENTER 

< 


_dtic 

J^HLECTE 


APR  0  3 1986 


D 


^PROrPnilir^^.Pf.^^SPECTION  SAMPLING 


BV  DEMING 


procedures  discussed 

by 


ORC  86-2 


January  1955 


UNIVERSITY  OF 


This  document  has  been  approved 
lot  public  release  and  sale;  its 
distribution  is  unlimited. 


CALIFORNIA 


BERKELEY 


4  2  048 


r, 


BAYESIAN  ANALYSIS  OF  INSPECTION  SAMPLING 
PROCEDURES  DISCUSSED  BY  DEMING 

by 

Richard  E.  Barlow  and  Xiang  Zhang 


ORC  96-2 


January  1986 


‘Department  of  Industrial  Engineering  nnd  Operations  Research,  University  of 
California,  Berkeley,  California. 


This  research  was  supported  by  the  U.  S.  Army  Research  Office,  under  Contract 
DAAG29-85*K-0208  with  the  University  of  Cniifornia.  Reproduction  in  whole  or  in 
part  Is  permitted  for  any  purpose  of  the  United  States  Government. 


MASTER  COPY 


SECURITV  CLMSiriCATION  OF  THIS  RACC  (VAan  Data  Bn>araA9 


REPORT  DOCUMENTATION  PAGE 


ORC  86-2 


A.  title  fanA  Sukiirfaj 

Bayesian  Analysis  of  Inspection  Sampling 
Procedures  Discussed  by  Deming 


T.  AUTHORfaJ 

Richard  E.  Barlow  and  Xiang  Zhang 


FOR  REPRODUCTION  PURPOSES  _  ^ 

M  h  Hi 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


S.  RECIHIENT’S  CATALOe  NUMBER 


B.  TYRE  OF  RERORT  A  RERIOO  COVEREO 

TECHNICAL  Report 


•  ■  RERFORMING  ORC.  RERORT  NUMBER 

ORC  86-2 


B.  contract  or  ORANT  NUMBERfa) 

DAAG29-85-K-0208 


B.  RERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Operations  Research  Center 

3115  Etcheverry  Hall 

University  of  California,  Berkeley 


RERORT  DATE 

January  1986 


IS.  number  OF  rages 

27 _ 


I  CaniraflAiB  Ofnca>  <S-  SECURITV  CL  ASS.  faf  M<a  noaft; 

Unclassified 

ISa.  DECLASSIFlCATION/DONNGRADING 


G.  mSTRIBUTION  STATEMENT  (ml  BUa  RaRarO 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ml  Cha  aSalracI  anfaraB  In  SfacA  SO.  If  Mtfmrml  fimm  ftmrmtt) 


IS.  surrlbmentarv  notes 


The  view,  opinions,  and/or  findings  contained  in  this  report  are 
those  of  the  authorCs)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy,  or  decision,  unless  so 


BTBAIkiiTUlSVzliaVSIf 


>  KBV  vmot  aid  Idwtfliy  ir  Mac*  maiDvr; 

Inspection  Sampling  Quality  Assurance 

Bayes  Decision  Analysis 

Influence  Diagrams 


SG  abstract 


’  mi  Maatlir  Sv  MacA  miSap 


(SEE  ABSTRACT) 


00  iSSTt.  un  COITION  OP  I  MOVES  IS  OOSOLCTE 


UNCLASSIFIED 
security  classification  c 


, 

I  /;■ 


BAYESIAN  ANALYSIS  OF  INSPECTION  SAMPLING  PROCEDURES 
DISCUSSED  BY  DENING^ 
by 

Richard  E.  Barlow  and  Xiang  Zhang 


ABSTRACT 

In  Chapter  13  of  Quality.  Productivity,  and 
Competitive  Position.  Doming  diecussee  in  detail 


inspection  sampling  relative  to  two  different  cost 
setups.  Because  of  their  practical  importance,  one  of 
these  models  was  examined  in  detail  from  a  Bayesian 

A  very 


point  of  view  by  Barlow  and  Zhang  (198^^. 
general  Bayesian  treatment  of  these  problems  is  provided 
by  Lorenzen  (1985|T  A  computer  program  is  also  provided 
by  Lorenzen  for  the  Bernoulli  case  as  well  as  inspection 
sampling  plans  which  reject  the  lot  when  the  number  of 
defective  items  is  sufficiently  large.  However,  as  we 
shall  see.  simple  rules  of  this  form  are  not  optimal  for 
all  of  Deming’s  models.  In  particular.  Deming’s 
discussion  of  inspection  rules  when  the  finished 
assembly  cannot  be  repaired  by  replacing  defective  units  7^ 

T 

is  of  this  type.  Deming  describes  this  model  as  value  II 
added  to  substrate'*.  For  example,  if  the  unit  is  a  bag^ 
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of  cement  powder,  then  it  cannot  be  recovered  from  the 
assembly  after  water  has  been  added. 

The  purpose  of  this  paper  is  to  analyze  and 
contrast  the  two  models  considered  by  Deming  from  a 
Bayesian  decision  analysis  point  of  view.  We  will 
analyze  the  models  using  influence  diagrams#[Howard  and 
Natheson  (1981).  Shachter  ( 1984) ] is"*a  relatively 
new  technique  for  studying  statistical  problems.  It 
will  be  useful  for  studying  inspection  sampling 
problems . 

2.  USING  INFLUENCE  DIAGRAMS 

We  will  explain  the  use  of  influence  diagrams  in 
terms  of  a  very  simple  model.  Suppose  we  are 
considering  whether  or  not  to  inspect  units  in  a  unique 
lot  of  size  N  of  such  units.  That  is.  there  are  no 
other  units  than  those  in  this  lot.  The  question  is 
whether  or  not  to  inspect  units  before  installing  them 
in  an  assembly.  Inspection  costs  ki  dollars  per  unit. 
If  a  unit  is  not  inspected  but  is  found  later  to  be 
defective  in  the  assembly,  the  cost  is  ka  dollars. 

The  decision  variables  are  n  .  the  inspection  sample 
size,  and  the  inspection  decision  after  all  n  units 
have  been  inspected.  Inspection  decision  d^  is  not  to 
continue  inspection.  Inspection  decision  di  is  to 
inspect  the  remaining  N-n.  Figure  2.1  describes  the 
problem  setup. 
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Figure  2.2  is  an  influence  diagram  related  to  our 


problem.  In  the  diagram,  circles  denote  unknown 


quantities  while  rectangles  denote  decision  nodes.  A 


double  circle  denotes  a  deterministic  or  logical  node. 


In  the  figure,  the  node  corresponding  to  Xa  .  the 


number  of  defectives  in  the  remainder  N-n,  is 


deterministic  since 


Xa  s  6  -  Xi  given  0  and  Xj  . 


Arrows  indicate  logical  and  statistical  dependence. 


Associated  with-  each  circle  is  a  conditional  probability 


distribution.  If  units  in  the  lot  are  Judged 


exchangeable  a  priori,  then  the  distribution  of  x^ 


given  0,  n.  and  N  is  hypergeometric .  An  initial  or 


prior  distribution  must  be  assigned  to  0  .  the  unknown 


number  of  defectives  in  the  lot.  Given  n  and  Xa  the 


total  cost  is  k^n  -t-  kaXa  if  the  decision  is  d^  .  If 


the  decision  is  d^  .  the  total  cost  is  k^N 


low  and  Zhang  (1985){>r  a  re.la.4e4  w4el. 


A  mathematically  convenient  distribution  model  for 


the  initial  distribution  of  0  is  the  Beta-Binomial.  0 


is  Beta-Binomial  if  0  given  p  is  Binomial (p.N)  and 


p  is  Beta(A,B)  where 

Beta(A.B)  JS  p^"'(l-p)®"' 


2WJ 


INSPECTION  OF  A  UNIQUE  LOT 
FIGURE  2.2 
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and  A.B  >  0  .  In  this  case  [cf.  Basu  and  Pereira 
(1982)] 

E[0-Xi !xi ,n,N]  =  (N-n) (A+Xi )/[A+B+n] 
and  the  decision  rule  di  is  optimal  if 
(A+Xi)/[A+B+n]  >  ki/ka  . 

Since  the  left  hand  side  is  increasing  in  Xi  and 
decreasing  in  n  <  N  .  we  have  that  nsN  is  optimal  if 
and  only  if 

A/[A+B+N-1]  >  kj/ka  . 

Unfortunately  the  Beta-Binomial  is  not  an 
appropriate  prior  distribution  if  we  strongly  believe  a 
priori  that  0  is  approximately  0^  (0  <  0^  <  N)  so 
that  the  coefficient  of  variation 

1/2 

a(0)/p(0)  <  C{l-V”)^®o^ 

For  the  Beta-Binomial 

«^(0)  =J[NA/(A+B)][1-A/{A+B)][(A+B+N)/(A+B+1)] 
and  p(0)  -  NA/(A+B)  s  0^  so  that  the  coefficient  of 
variation  satisfies 


a(0)/4(0) 


(l-0^/N)  (A+B+N) 

0^  (A+B+1) 


>  J  (l-0/H)/0^ 


In  this  case  d  may  be  best  if  x.  is  sufficiently 
o 

large  so  that  the  decision  rule  is  reversed  from  the 


Beta-Binomial  case. 

In  the  appendix  the  solution  strategy  for  the 
problem  with  arbitrary  prior  is  described  in  terms  of 
arc  reversals  and  node  eliminations.  To  determine  the 
sample  size,  n  . 

Minimize  I  Min{k}n  ksECxs  Sxi  .n.N) .  kiN}  p(Xi!n.N)  . 
n  Xj 

Knowing  the  decision  rule  (as  in  the  Beta-Binomial  case) 
simplifies  the  computation  by  allowing  the  first  minimum 
computation  before  the  summation. 

3.  DENING'S  INSPECTION  SAMPLING  MODELS 

Units  from  a  lot  of  size  N  are  Installed  in  an 
assembly  at  a  point  in  the  production  line.  The 
question  is  whether  or  not  to  inspect  the  unit  (at  cost 
$ki/unit)  before  Installing  the  unit  in  an  assembly.  If 
the  unit  is  not  inspected  but  is  defective,  then  this 
will  be  discovered  at  assembly  test.  A  cost  $k2  >  $kt 
will  then  be  incurred.  In  model  1.  the  defective  unit 
found  at  assembly  test  will  be  replaced  by  a  good  unit. 
In  model  2,  this  is  not  possible  and  the  assembly  will 
be  sold  at  a  reduced  price.  In  this  case  ka  is  the 
loss  Incurred.  In  Deming’s  models  all  lots  are  Judged 
exchangeable  a  priori  so  that  a  binomial (p.n) 
probability  function  is  used  for  the  number  of  defective 
units  in  the  lot.  The  unknown  percent  defective  is 
denoted  by  p. 


In  both  models  we  consider  inspecting  n  N  units. 
Let  Xi  be  the  number  defective  in  the  sample  of  size 
n.  Since  known  defective  units  will  not  be  Installed, 
an  additional  units  are  inspected  until  good  units 

are  found  and  installed.  Since  the  vendor  will  replace 
defective  units  at  the  vendor’s  cost  we  may  consider  the 
yi  units  that  are  inspected  to  come  from  a  separate  and 
very  large  lot  provided  by  the  vendor.  The  sampling 
inspection  cost  is  ki^n-fy^]  . 

After  inspecting  n  units  we  consider  only  two 
possible  decisions  for  both  models,  namely 
d^:  STOP  inspecting  units 

df!  INSPECT  the  remaining  N-n  . 

In  model  1.  the  defective  units  found  in  the 

remaining  N-n  units  at  assembly  test  are  replaced  by 
good  units.  To  do  this  an  additional  units  need  to 

be  inspected  to  find  good  units.  The  y^  units 

Inspected  also  come  from  the  separate  lot  provided  by 
the  vendor.  The  cost  of  defective  units  found  at 
assembly  test  is 

kaXa  *  kiya  . 

The  logical  and  statistical  dependency  between  the 
unknown  quantities  and  decisions  are  shown  in  Figure 
3.1.  The  random  quantity,  p  ,  is  the  unknown  percent 
defective  in  many  lots,  all  Judged  exchangeable. 
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Figure  3.1  makes  explicit  the  well  known  fact  that 
Xi  and  X2  given  n  and  p  are  statistically 
independent.  Table  3.1  describes  the  variables  used  in 
Figure  3.1. 
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TABLE  3.1 
MODEL  1 


.n) 


di(xj ,n) 


Sample  Size 

Decision  to  STOP  inspection  given 
(xi .n) 

Decision  to  INSPECT  remaining  N-n 
given  (xi.n) 


DECISION 

DESCRIPTION  \ 

VARIBLES 

u 

FIGURE  3.1 
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Model  1  was  examined  in  detail  in  Barlow  and  Zhang 
(1985).  Let  Zi  be  an  indicator  for  the  i-th  unit 
which  is  1  if  the  unit  is  defective  and  0  otherwise. 
Given  n  and  ,  d^  is  optimal  if 

E(£(Zi ip) !xi .n)  i  ki/ka  .  The  assumption  made  was  that 
only  Xi  and  n  were  known  at  the  time  of  the 
inspection  decision.  Since  yi  .  the  additional 
inspections  required  is  also  presumably  known,  there  is 
a  loss  of  Information.  If  both  n  and  y^  are  known 
at  the  time  of  the  inspection  decision  then  the  critical 
number.  c(n)  .  may  exceed  n  in  this  case.  Other 
than  this,  all  the  results  in  Barlow  and  Zhang  (1985) 
carry  over  to  the  case  that  yi  is  known  at  the  time  of 
the  inspection  decision. 

4.  DENING’S  MODEL  FOR  VALUE  ADDED  TO  SUBSTRATE 

Work  is  done  on  incoming  material,  the  substrate. 
The  finished  product  will  be  classed  as  first  grade  or 
second  grade  or  third  grade  or  scrap.  Let  ka  be  the 
average  loss  from  downgrading  the  final  product  or  for 
scraping  finished  items. 

If  the  incoming  unit  is  non-defective,  then  it  will 
result  in  a  final  product  of  first  grade;  otherwise,  if 
it  is  defective,  and  if  it  goes  into  the  production 
line,  then  a  final  product  of  downgrade  or  scrap  will  be 
produced.  Every  final  product  is  subject  to  assembly 
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grading.  Figure  4.1  is  an  influence  diagram 

representation  of  our  problem.  In  this  case  there  is  an 

arrow  from  inspection  decision  to  y2  since  yg  =  0  if 

decision  d  is  taken, 
o 

The  optimal  inspection  decision  rule  can  be 

calculated  for  every  sample  size  n  in  terms  of  the  costs 

ki .  ks  and  the  prior  distribution  for  p. 

Lemma  1 .  Suppose  n  units  are  inspected  and  Xi  are 

found  defective.  The  optimal  inspection  decision  rule 

in  this  case  is  d^  if 

E[Var(Zj!p)  !  Xj.n-1]  i  k^/ka 

where  expectation  is  with  respect  to  the  posterior 

distribution  of  p  given  Xf  and  n-1  . 

Since  p(l-p)  i  1/4  for  all  0  ^  p  ^  1  ,  it 

follows  from  Lemma  1  that  if  ki/ka  >  1/4  .  then  d  is 

o 

optimal  for  all  possible  0  ^  Xj  ^  n  ^  N  ;  i.e.  n  s  0 
is  optimal  in  this  case. 

For  a  Beta(A.B)  prior  and  given  Xt  and  n  .  di 
is  optimal  if 

(A+Xi  ){B+n-l-Xi  )/[(A-»-B+n)(A+B+n-l)]  i  ki/ka  • 

Hence,  given  x^  and  n  .  di  is  optimal  if 

X,  e  [cj(n)  .  ca(n)] 
where  Ci(n)  and  Ca(n)  are  given  by 


L5 


CaCn) 


-b  ♦ 

2a 


4ac) 


where  a  =  1,  b  =  -(B+n-l-A), 

c  =  -A(B+n-l)  (A+B+n)(A+B+n-l  )ki/k2  • 

The  intuitive  reason  for  this  rule  is  as  follows. 

If  the  expected  variance  is  snail,  then  we  expect  that 
p  is  either  close  to  zero  or  close  to  one.  If  we 
believe  p  is  close  to  0  so  that  there  will  be  few 
defectives,  then  there  is  no  point  to  inspection.  On 
the  other  hand,  if  we  believe  p  is  close  to  one,  then 
there  will  be  so  many  defective  units  that  the  cost  of 
inspecting  units  to  find  a  good  unit  to  replace  a 
defective  unit  before  asseably  will  tend  to  exceed  the 
cost  due  to  downgrading  the  final  assembly  even  if  it  is 
defective . 

This  decision  rule  is  in  contrast  to  the  decision 
rule  for  Model  1  where  we  inspect  all  units  if  the 
posterior  expected  probability  that  a  unit  is  defective 
is  sufficiently  large. 

Proof  of  Lemma  1.  Suppose  we  inspect  n  units  and  x^ 

are  found  defective.  The  expected  cost  given  n  and 

Xi  if  we  take  decision  d^  is 

1 

c(d-)  »  (N-n)  J  k2Pv(p!X|.n}  dp  . 

®  0 

If  decision  d^  is  taken,  then  the  expected  cost  is 

1  .1 

C(d»)  «  (N-n)  jr  ki(l-p)  ‘w(p!x,.n)  dp  . 

0 
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Consider  the  difference 


C{d  )  -  C(d,)  =  (N-n)k2  /  [p  -  k^/ka  ( l-p)“' ]ir(p  txj  .  n)  dp 
°  0 

1 

=  (N-n)k2  /  (1-p)  [p(l-p)  -  ki/k2]»(p!xi .n)  dp  . 

0 

Since  ir(p!xi.n)  is  proportional  to  P**  ( l“P)”~*‘x(p) 

1  . 

C(d  )  <  C(d.)  iff  /[p(l-p)  -  k,/k2]p*‘(l-p)"“^"*‘x(p)  dp  <  0 
®  0 

1 

iff  X  p(1“p)t(p!xi ,n-l)  dp  <  k^/ka 
0 

1 


iff  X  Var(Zj :p)»(p:x*.n-l)  dp  <  kj/ka 
0 

iff  E[Var(Zi!p)!xt.n-l]  <  kj/ka  .  QED 


Lemma  2.  E[Var (Zj !p) Sxi .n-1]  is  unimodal  in  Xi  (n) 
when  n  (xi)  is  fixed. 

1 

Proof.  E[Var(Zt  !p)  !x,  .n-1]  «  X  p(  l-p)ir(p  Sx^  ,n-l )  dp  . 

0 

Note  that  ir(p!xi.n-l)  can  be  written  as 
g(xi  ,n) ( l-p)p*‘ ( 1-p)"  .  It  is  well  known  that 

p^*(l-p)"~^*  is  a  totally  positive  function  of  Xi  and 
p  and  reverse  regular  function  of  n  and  p  [Karlin 
Theorem  2.1,  page  18  (1968)].  It  follows  that 


|T 
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T(p;xi,n-1)  is  totally  positive  in  Xi  and  p  and 

reverse  regular  in  n  and  p  since  g(xi.n)  >  0  . 

Using  the  sign  variation  diminishing  properties  of 

sign-regular  functions  [Section  3,  Chapter  5.  Karlin 

(1968)3  we  can  finish  the  proof.  To  do  this  consider 

1 

h(xi.n)  -as/  [p(l-p)  -  a]w(p!xt .n-1)  dp 
0 

1  2 

■  X  [-P  ♦P-»3»(p!*i .n-1)  dp 
0 

for  a  fixed  real  a  .  The  quadratic  function  can  have 
at  most  two  sign  changes  in  [0.1]  •  The  sign  must 
change  from  negative  to  positive  to  negative  if  there 
are  exactly  two  sign  changes.  It  follows  that 
h(xi.n)  -  a  as  a  function  of  Xi  (n)  when  n  (xi)  is 
fixed,  changes  sign  at  most  two  times  when  Xi  (n)  runs 
from  0  (Xi)  to  n  (N)  .  Furthermore  h(xi.n)  as  a 
function  of  Xj  changes  sign  in  the  same  order  if  there 
are  two  sign  changes.  As  a  function  of  n  it  changes 
sign  in  the  opposite  order  if  there  are  two  sign 
changes.  Since  a  is  an  arbitrary  real  number,  this 
implies  the  unimodality  properties  as  stated.  QED 

The  figures  below  illustrate  h(xi.n)  as  a 
function  of  Xi  for  fixed  n  and  also  as  a  function  of 
n  for  fixed  Xi  . 
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Theorem.  For  Model  2  and  given  n  and  the  optimal 

inspection  decision  is  d^  if  and  only  if 

Xi  e  [C4(n).  ca(n)]  . 

where  Ci(n)  and  C2(n)  are  solutions  to 
1 

/  p(l-p)»{p!c.n-l)  dp  s  0  . 

0 

This  is  true  for  any  prior  distribution  for  p. 

The  proof  is  an  immediate  consequence  of  Lemmas  1 
and  2 . 

5.  CONCLUSIONS 

Influence  diagrams  are  a  useful  technique  for 
displaying  the  logical  and  statistical  dependence 
between  unknown  quantities  and  decision  variables  in 
Inspection  sampling  problems.  They  can  also  be  used  to 
develop  a  computational  solution  strategy.  The  solution 
strategy  is  simplified  if  the  form  of  the  optimal 
inspection  decision  rule  can  be  calculated.  However, 
the  "obvious”  inspection  rule  of  inspecting  all 
remainders  in  a  lot  when  the  number  of  defectives  in  a 
sample  is  sufficiently  large  is  not  always  correct. 

This  has  been  shown  for  two  cases;  namely  the  case  of  a 
unique  lot  when  we  have  a  strong  prior  for  0  (i.e.  the 

prior  coefficient  of  variation  is  sufficiently  small) 
and  for  the  Doming  Model  2  corresponding  to  "value  added 
to  substrate.” 
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